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ABSTRACT 



Aims. We present a study of the neutron star X-Ray Transient Cen X-4. Our aim is to look for any evidence of irradiation of the 
companion with a detailed analysis of its radial velocity curve, relative contribution of the donor star and Doppler tomography of the 
main emission lines. 

Methods. To improve our study all our data are compared with a set of simulations that consider different physical parameters of the 
system, like the disc aperture angle and the mass ratio. 

Results. We conclude that neither the radial velocity curve nor the orbital variation of the relative donor's contribution to the total 
flux are affected by irradiation. On the other hand, we do see emission from the donor star at Ho- and Hel ,45876 which we tentatively 
attribute to irradiation effects. In particular, the Ha emission from the companion is clearly asymmetric and we suggest is produced 
by irradiation from the hot-spot. Finally, from the velocity of the Hel /15876 spot we constrain the disc opening angle to a = 7° - 14°. 
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1. Introduction 

Soft X-Ray Transients (SXRTs), a subclass of Low Mass X-Ray 
Binaries, are interacting binary systems which alternate short 
periods (weeks to months) of high X-ray luminosity with long 
(several years) intervals of quiescence. During quiescence, such 
transient systems, are very faint in X-rays (10 32 - 10 33 erg s" 1 ) 
and their optical luminosity drops by as much as 6-7 mag, giv- 
ing a unique opportunity for the study of companion stars (see 
Campana et al. 1998 for a review). 

Cen X-4 is a well-known neutron star X-ray tran- 
sient, from which two bright X-ray outburst were de- 
tected in 1969 HConner, Ev ans & Belian 1969) and in 1979 
( |Kaluzienski, Holt & Swank 1980). From then on it remained 
in the quiescent state. With an evolved 0.3 M Q , K3-K7 star 
(Shahbaz, Nay lor & Charles 1993|ITorres et al. 20 02 D'Avanzo 
et al. 2005), Cen X-4 is one of the brightest SXRTs in 
quiescence, with V = 18.7 and low interstellar absorp- 
tion (Ay = 0.3 mag). A 15.1 hr orbital period was de- 
termined from the sinusoidal variation of the optical light 
curve (Cowley et al. 1988; Chevalier et al. 1989; McClintock 
& Remillard 1990). Both the optical spectrum and emission 
line maps obtained w ith the Doppler tomography technique 
( Marsh & Home 1988) show clear evidence of disc emission 
JTorres et al. 20021 IrTAvanzo et al. 20051 ). 

In general, the effect of heating of the secondary star leads to 
a change in absorption lines strength, with the consequent effect 
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to move the effective centre of the secondary away from the cen- 
tre of mass of the star. This causes phase-dependent distortions 
of the radial velocity from a sinusoidal fit, leading to incorrect 
determinations of the radial velocity amplitude. 

In our previous work on Cen X-4 evidence of irradiation was 
proposed in light of a residual emission of Ha and Hel from 
the companion star. We present here a search for evidence of 
irradiation in Cen X-4, comparing our observations presented 
in D'Avanzo et al. (2005) with the simulations obtained using 
the code presented in Mufioz-Darias, Casares & Martinez-Pais 
(2005). 

2. The radial velocity curve of Cen X-4 

In D'Avanzo et al. (2005) we presented the results of the first 
extensive and contiguous orbital coverage of Cen X-4. In Fig.^ 
the radial velocity curve we obtained for Cen X-4 is shown, to- 
gether with the residuals to the sinusoidal fit. To overcome any 
spurious effect, we have carefully corrected our spectra for any 
telescope flexure, using telluric lines. It is clear from Fig.^that 
there are deviations from the best fit; such deviations are of the 
order of a few km s _1 . 

In order to understand if this behaviour could be considered 
as an indication of irradiation we have compared our observa- 
tional result with a simulation obtained with the code developed 
by Munoz-Darias, Casares & Martinez-Pais (2005). In this simu- 
lation we model the effect of irradiation in radial velocity curves 
derived from the Doppler shift of absorption lines formed by the 
companion. These features are quenched in the heated face of 
the donor and hence the observed Zf-velocity provides only an 
upper limit to the real K%. This deviation roughly depends on 
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Fig. 1 . Radial velocit y curve for Cen X-4 
(from lD 7 Avanzo et al. 20 05 1 obtained through cross-correlation 
of absorption lines in the range 5930 - 6750A with late-type 
template stars. Two phases are shown for clarity. The best 
sine-wave fit is also shown. Residuals from the best fit are 
shown in the lower panel. 

the X-ray luminosity, the orbital separation and the opening disc 
angle (a) which partly obscures the companion. Here, we have 
only considered the limit case in which the absorption lines are 
totally quenched on the irradiated regions of the donor star which 
are not shadowed by the disc. 

We have produced two curves (both showed in Fig.|2ji, repre- 
senting the irradiated and the non-irradiated case. Residuals are 
shown at the bottom of the figure. It is immediately clear that the 
main effect of irradiation from the neutron star is the rise in the 
velocity semi-amplitude, as shown. Such a behaviour is clearly 
observed in the radial velocity curves of systems where the sec- 
ondary st ars are exposed to irradiation from the c ompact object 
(like, e.g., IFriendetal. 1990allFriend et al. 1990bl 

In light of the results of our simulations, we can conclude 
that the radial velocity curve does not show any evidence of ir- 
radiation. The expected rise in the velocity between phases 0.25 
and 0.75 (i.e. when the observer sees the inner face of the com- 
panion) is not present. Although deviations from the sinusoidal 
fit seem to be present in our observed curve these are signifi- 
cantly smaller (given the error bars) than predicted by the irra- 
diation model, and most important, these deviations are in the 
opposite direction to what is expected in an irradiation scenario. 

3. The optical light curve of the companion 

If the inner side of the secondary star is exposed to irradiation 
this would also affect the modulation of the contribution of the 
companion to the optical luminosity. In D'Avanzo et al. (2005) 
we determined the variation of the fractional contribution of light 
from the secondary (named factor/) for 30 phase bins. The result 
was that/ is not constant at all the orbital phases, but modulated 
with two unequal minima at phase (i.e. at superior conjuction) 
and 0.5 and two nearly equal maxima at phases 0.25 and 0.75 
(i.e. when the observer sees the sides of the companion with the 
maximum surface projected area). This behaviour is reminiscent 



Fig. 2. Simulated radial velocity curves for Cen X-4 for the ir- 
radiated (stars) and non irradiated (triangles) case. Differences 
between the two models are shown in the lower panel. 
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Fig. 3. Variations of the fractional contribution of the companion 
compared with synthetic models for a = 0° (dotted line), 7° 
(solid line) and 14° (dashed line). 



of the classical ellipsoidal modulation (see e.g. McClintock & 
Remillard 1990; Shahbaz et al. 1993) caused by the changing 
visibility of the tidally distorted companion as seen by the ob- 
server. However, we note the unusually large difference in depth 
between phases and 0.5, requiring an additional effect like, 
e.g. the heating of the inner hemisphere of the companion by 
X-ray irradiation from the neutron star. In order to investigate 
this scenario we compare the variation of the / factor (coad- 
ded into 14 phase bins to increase the S/N ratio) with a set of 
synthetic curves obtained with the code developed by Munoz- 
Darias, Casares & Martinez-Pais (2005). This is presented in 
Fig.fJ 
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These simulations show the absorption line flux as function 
of the orbital phase for three different flaring angles (a). It is 
well established that the temperature of the outer layers of the 
companion can be upset by external heating (e. g. Brett & Smith 
1993), which removes the vertical temperature gradient of the 
heated star and quench the absorption lines within the irradiated 
regions. In order to compute our curves we have considered the 
limit case in which absorption lines are totally quenched (i. e. 
not present) on the irradiated regions of the companion (see also 
Shahbaz et al. 2000). As it is clear in the plot we cannot reconcile 
the models and the data for any value of a. If the opening angle 
is low, the absorption line is quenched upon the whole irradiated 
hemisphere of the companion and hence the absorption line flux 
is very low at orbital phases close to 0.5 (dotted line). However, 
we can not reproduce the lower / factor found for orbital phases 
around 0. On the other hand, for a > 14° (dashed line) the ac- 
cretion disc completely screens out the companion star, the effect 
of irradiation disappears and we obtain the ellipsoidal modula- 
tion associated with the changing visibility of the tidally dis- 
torted companion star. If we consider the X-ray luminosity of 
Lx ~ 10 32 erg s _1 reported in Campana et al. (2004) and the 
orbital parameters of the system, we estimate using our binary 
code that the effective temperature of the irradiated regions of the 
donor is only - 30 K higher than in the non-irradiated regions. 
Hence, the donor's atmosphere should almost not be affected by 
heating and we expect an / factor behaviour close to the latter 
case in which irradiation effects are not present. Therefore, as 
result of our simulations, the irradiation from the inner regions 
of the accretion disc does not appear to be responsible for the 
orbital behaviour of the / factor. 



4. Ha and Hel emission from the companion 

In D'Avanzo et al. (2005) we presented Doppler maps of Cen 
X-4 showing clear evidence of residual Ha and Hel ,15876 
emission associated with the companion star. As opposed to 
Hel ,15876, the Ha spot is clearly offset with respect to the 
vertical V y axis. As reported in D'Avanzo et al. (2005) we have 
computed the centroids of these spots obtaining: 

Ha: (V x ,Vy) = (34.1 + 11.8,97.7 ± 11.8) km s" 1 

Hel ,15876: (V x ,V y ) = (10.3 ±11.8, 122.8 ± 11.8) km s" 1 

This clearly demonstrates that the offset of the Ha spot is signif- 
icant at ~ 3 - cr and it leads phase by 0.05 + 0.03 cycles. The 
S-wave component responsible for the Ha spot is clearly seen in 
the individual spectra. Therefore, in order to reduce the uncer- 
tainty in the position of the region of the companion emitting Ha 
photons we have performed a multi-Gaussian fit to the Ha pro- 
files, coadded into 30 phase bins. Such operation enabled us to 
isolate the Ha component emitted from the companion for each 
phase bin, and to determine with good precision its centroid. 
As a result we obtained a sinusoid-like velocity curve, shown 
in Fig. whose semi-amplitude (equal to 91.1 + 3.0 km s _1 ) 
is in agreement with the value of Vy obtained from the Doppler 
map, but it is more precise and gives the velocity of the Ha emis- 
sion region of the companion. Moreover, we also note that Fig.@] 
shows also an indication of a phase shift (equal to 0.02 + 0.01) 
which is in agreement to the shift in the V x direction of the Ha 
spot in the Doppler map. 

Unfortunately the lower S/N ratio of the Hel ,15876 line does 
not enable us to do an accurate multi-Gaussian fit to isolate the 
companion component, so we used the position given from the 
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Fig. 4. Velocity of the Ha emission from the companion at dif- 
ferent orbital phases. Two cycles are shown for clarity. 

Doppler map. In this case no phase/velocity shift is present in 
the x direction. 

In D'Avanzo et al. (2005) we measured an equiva- 
lent width for the Ha residual of 4.4 + 0.5 A which 
could not be explained only invoking chromospheric ac- 
tivity from the companion. Hence, we proposed that the 
Ha spot may be triggere d by irradiation from the com- 
pact star (for details see ID'Avanzo et al. 2005> . This sce- 
nario predicts lower excitation lines to be formed at higher 
Roche lobe latitudes because of variable photoelectric ab- 
sorption by the outer disc rim i.e. lower energy photons are 
more efficiently absorbed and hence lower excitation lines 
arise from regions close to the poles. This has been previ- 
ously seen in the cataclysmic variable IP Peg (Harlaftis 1999, 
Morales-Rueda, Marsh & Billington 2000 1 and more recently in 
U Gem ( |Unda-Sanzana, Marsh, & Morales-Rueda 20 06 1. This 
scenario, however, does not agree with our observed veloc- 
ities for Ha and Hel emission. Moreover, with a velocity 
semi-amplitude of 91.1 + 3.0 km s , the Ha residual emission 
arises from the vicinity of the Lj point which must be efficiently 
shadowed from the UV photons by the disc. This, together with 
the asymmetric distribution of the Ha spot over the Roche lobe, 
points to the hot-spot as a possible source of the Ha photoionis- 
ing radiation. 

5. Discussion 

The discrepancies between the observed and the simulated ra- 
dial velocity curve enable us to conclude that our radial velocity 
curve of Cen X-4 does not contain any evidence of the pres- 
ence of an irradiated companion and that the measured semi- 
amplitude of the curve must correspond to the "real" K2 velocity 
of the secondary star. This is also confirmed by new high quality 
UVES data at 7 km s _1 resolution (Casares et al. 2006, in prepa- 
ration). Therefore, no "K-correction" seems to be necessary for 
this system. 

The interpretation of the orbital modulation of the optical lu- 
minosity contribution of the companion is more trivial. We note 
here that the /-factor yields not the lightcurve of the companion 



4 



P. D'Avanzo et al.: Irradiation in Cen X-4 (RN) 



E 



1 40 




,,,,,,,,,,, 


. 


130 




/ / 




120 




q=0.17X JL^ 11 


J 


110 
















1 00 
90 









2 4 6 8 10 12 14 

DISC FLARING ANGLE ( s ) 

Fig. 5. Determination of the disc flaring angle in Cen X-4. The 
solid lines represent the result of our model for the considered 
mass ratio values (q ~ 0.17 and q ~ 0.23). The dashed lines 
represent the value K em = 122.8 ±11.8 km s measured from or 
Doppler maps. The possible disc flaring angle range is marked 
with a solid line. 



star but the relative contribution to the total flux. Therefore, in- 
stead of neutron star irradiation, a more plausible explanation for 
the deep phase 0.5 minimum would simply be an increase in the 
disc brightness around this particular phase i.e. a non-uniform 
emissivity distribution which might trace azimutal variations of 
the disc density. 

Neutron star irradiation seems to be ruled out also as respon- 
sible for the Ho- emission from the companion. In fact, as we 
have seen in the previous section, the Ha residual arises from a 
region of the companion which should be affected by the shadow 
of the disc. 

In light of this, hot-spot irradiation seems to be a natu- 
ral alternative explanation for the intrinsic Her emission from 
the companion, because in such case, no shielding is to be 
considered. In fact, as we see from our Doppler maps, bright 
anisotropies are present on the accretion disc. We propose that 
the Ha emission concentrates in the trailing side of the Roche 
lobe because the leading side is shielded by gas stream. Hot- 
spot heating has been previously invoked to explain asymmetric 
emission patterns observed in several cataclysmic variables (e.g. 
in IX Vel, see Beuermann & Thomas 1990). Furthermore, evi- 
dence for a hot-spot or splash point in Cen X-4 was presented 
by McClintock & Remillard (1990) whereas our Hel Doppler 
images (Fig. 4 in D'Avanzo et al. 2005) also show evidence of 
bright spots coincident with the interaction between the outer 
disc rim and the gas stream. Alternatively, asymmetric heating 
may be amplified by the action of circulation currents induced 
through Coriolis forces (Davey & Smith 1992, Martin & Davey 
1995). A direct consequence of hot-spot (if present) or, in gen- 
eral, anisotropy irradiation is the shift, in the V x direction, of the 
Ha emission region on the companion, exactly as observed. 

On the contrary, the higher V y and the absence of any 
phase/velocity shift in the x direction, should suggest that the 
Hel emission region on the companion lies in a region of the 
Roche Lobe that is out from the shadow of the disc and directly 



exposed to illumination by the neutron star. We can use this hy- 
pothesis to put some strong constraints on the aperture angle a of 
the accretion disc. To do this, we have used the K-correction val- 
ues tabulated in Munoz-Darias, Casares & Martinez-Pais (2005). 
In this work, the relation K em /K2 is computed as function of the 
mass ratio (q — j±) for different values of a. Where Ki is the 
radial velocity associated with the mass center of the companion 
and K em is the velocity of an emission line formed on the heated 
face of the donor. Since q and K2 are known for Cen X-4 we 
have used the observed Kn e i emission line velocity to constrain 
the value of a. As is explained above, the low S/N ratio of this 
line does not allow us to isolate the companion component and 
hence we have used the K em = 122.8 +11.8 km s _I obtained 
from the Doppler map. In Fig. 5 we show that using this value of 
K em and considering 0.17 < q < 0.23 we directly obtain a > 7°. 
Moreover, if we take into account the maximum value of a to 
the companion be irradiated for this q, we can further constrain 
the opening disc angle between 7° < a < 14° (Fig.[5}- 

6. Conclusions 

We have examined the Ha and Hel /15876 bright spots associated 
to the companion star in Cen X-4. The Hel spot is well centered 
on the line joining the two stars and we use its velocity to con- 
strain the disc flaring angle to be 7-14°. On the other hand, the 
Hq, spot concentrates around the LI point and is shifted towards 
the leading side of the Roche lobe. We tentatively explain this 
asymmetry through irradiation from the hot-spot and shielding 
by the gas stream. 

We have also searched for evidence of irradiation effects in 
the radial velocity curve of the photospheric absorptions, by 
comparing our data with model simulations, but found none. 
This lead us to conclude that the donor's metallic lines are not af- 
fected by irradiation and hence the orbital parameters presented 
in previous papers (Torres et al. 2002, D'Avanzo et al. 2005) are 
safe from bias. 
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